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Qj ' Abstract 

^ : 

' We present numerical simulations of fragmentation of the Affleck-Dine conden- 

' sate in two spatial dimensions. We argue analytically that the final state should 

}—i ' consist of both Q-balls and anti-Q-balls in a state of maximum entropy, with 



most of the balls small and relativistic. Such a behaviour is found in simulations 
on a 100x100 lattice with cosmologically realistic parameter values. During frag- 
mentation process, we observe filament-like texture in the spatial distribution 
of charge. The total charge in Q-balls is found to be almost equal to the charge 
in anti-Q-balls and typically orders of magnitude larger than charge asymmetry. 
Analytical considerations indicate that, apart from geometrical factors, the re- 
sults of the simulated two dimensional case should apply also to the fully realistic 
three dimensional case. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) has several flat directions 
where the scalar potential is nearly identically zero During inflation, squark and 
slepton fields will fiuctuate freely along the flat directions, forming Affleck-Dine (AD) 
-condensates 0. The state of lowest energy, however, is not the AD-condensate but 
a non-topological soliton, the Q-ball 0, ^, which carries a non-zero baryonic (B-ball) 
or leptonic (L-ball) charge. The instability is induced by the spatial perturbations 
[^, which are naturally present in the condensate because of quantum fluctuations 
during inflation. The fragmentation process and the properties of the resulting Q-balls 
depend on SUSY breaking: if gauge mediated, Q-balls will be large and completely 
stable [^, 0] and form along every flat direction^ and could be detectable even today 
whereas if SUSY breaking is gravity mediated, Q-balls will be unstable but 
long-lived enough to have a host of interesting cosmological consequences [0, |n[]. In 
the gravity mediated case the formation of Q-balls is also a generic feature in all but 



few flat directions 12 



Fragmentation of the AD-condensate involves highly non-linear dynamics and is 



therefore quite complicated |[13|, The condensate lumps must somehow loose the 



extra energy when settling down to a Q-ball conflguration. This they may do by 
radiating out quanta of AD-scalars or smaller condensate lumps, but the lumps may 
also experience frequent collisions. Since for a flxed charge the energy of a Q-ball is 
also flxed, fragmentation dynamics and the flnal distribution of Q-balls will obviously 
depend on the initial energy and charge density of the condensate. In principle, this is 
a free parameter, although the observed baryon asymmetry can be argued to indicate 
a natural order of magnitude for the ratio. 

In Ref. the relaxation of a single spherical condensate lump was followed nu- 
merically, but sphericality is likely to be a far too constraining assumption. Moreover, 
the condensate fragments are not isolated but interact with each other. The general 



features of Q-ball collisions have been studied in Refs. ||T5|, |T^, |T7|; the cross section in 
the gravity mediated case has been shown to depend on the relative phases of the col- 
liding Q-balls 0]. Q-ball formation has been seen in 3-d lattice studies both in gauge 
mediated []TB| and gravity mediated cases but only for a single energy-to-charge 
ratio of the condensate in each case. As we will show in this paper, natural choices for 
the initial energy-to-charge ratio necessarily lead to a copious formation of not only 
Q-balls but also of anti-Q-balls, which carry a net negative charge. Such a conclu- 
sion can be reached already by relatively simple analytical arguments which assume 
that the condensate lump reaction rates are fast enough to drive them into a state of 
maximum entropy. From the resulting equilibrium distributions it then follows that 



^For another variant of gauge mediated Q-balls, see 



the number of Q-balls is almost equal to the number of anti-Q-balls so that the total 
positive (negative) charge in Q-balls is typically much larger than the actual charge 
asymmetry. Indeed, as we shall discuss in this paper, detailed numerical simulations 
show this behaviour. 

The paper is organized as follows. In Section 2 we derive the Q-ball and anti- 
Q-ball number density distributions by assuming that the condensate lumps actually 
equilibrate; this is shown to be a self-consistent assumption. We also show that as far 
as the distributions are concerned, there is no qualitative difference between two and 
three spatial dimensions. In Section 3 we present the results of numerical simulations 
of the condensate fragmentation in 2+1 dimensions. We follow the time evolution of 
the condensate field, show it breaking up into Q-balls and anti-Q-balls, and extract 
the number density distributions from the numerical data. In Section 4 we present 
our conclusions. 



2 Distributions of Q-balls and anti-Q-balls 

2.1 AD-condensate fragmentation 

In the cosmological context the formation of Q-balls via fragmentation of the AD- 
condensate has to be characterized by taking into account the expansion of the uni- 
verse. The random perturbations on the homogeneous AD-condensate field generated 
during de Sitter -phase of the inflation first grow linearly as long as they remain small 
compared to the AD-field. When these linear modes are large enough, they become 
non-linear mode by mode and their dynamics becomes much more complicated. During 
this non-linear era the actual Q-balls are formed as the system finds the energetically 
most favourable configuration ending up with a particular system of Q-balls and, in 
general, radiation. 

To be able to quantify the dynamics of the Q-balls one has to specify the potential. 
It consists of various terms appearing from the effective field theory of the particular 
case. To lowest order, a generic general potential reads 

VH{<l>) = V{<l>)-gH'\<!>\\ (1) 



where the second term, depending on the Hubble rate H = a/a, where a{t) is the scale 
factor, originates from higher-order operators in the Kahler potential (with typically 
g ^ 1), while V includes all terms of the flat direction. 

Of particular interest are flat-direction potentials of the form 

V{^) = m'lmi + log(^)) + (2) 



arising in the gravity mediated SUSY breaking scenario with a d=4 flat direction, and 

i 

V(»)=™*ln(l + ^)H-^H« (3) 

in the gauge mediated SUSY breaking scenario. The mass scale m is given by the 
SUSY breaking scale; typically m ~ (9(100) GeV. We have here omitted all terms 
which violate the quantum number Q. These are surely needed for generating the 
net charge of the AD-condensate, but the Q- violating terms are negligible at the time 
when the condensate finally begins to fragment. Therefore we may simply assume as 
an initial condition that at the onset of fragmentation the AD-condensate has some 
non-zero initial charge Q. 

The initial state of the AD-condensate is determined by two parameters, the energy 
and charge density of the condensate. For a fixed charge Q the critical parameter 
determining the Q-ball distribution is the ratio of energy density p to charge density 

One should note that according to our assumptions the charge density is a conserved 
quantity whereas the energy density is not. However, a simple calculation reveals that 
in a matter dominated universe the energy pumped into the scalar field by gravitation 
is very small compared to the initial energy of the condensate, so that in practice also 
p is conserved. Hence, the ratio x is a good parameter to describe the system. It is 
clear that the larger x is the more energy should be stored as radiation, kinetic energy 
of Q-balls, anti-Q-balls or some combination of these. 

To get an impression of the general features of the Q-ball formation process, it is 
useful to study instabilities in both gauge and gravity mediated scenarios. By writing 
the AD-condensate in terms of the modulus and the phase, 

$ = 0e^^ (5) 

one can formalize the requirements for the occurrence of the instabilities 5$ in the 
model. When the Q-balls begin to form, the instability band of fluctuations in the 
gauge mediated case is essentially given by |]18| 



* (6) 



a{t) 



where it is assumed that the field (p is large. This implies, that the size of the insta- 
bilities is bounded below by 

1 l$l 



In a cosmological context where the field is large this implies that the Q-balls will 
be very large. Indeed, this can be seen in the study by Kasuya and Kawasaki 
where they establish the formation of Q-balls using 3D lattice simulations. In their 
simulation one can see a single large Q-ball containing virtually all of the initial charge 
set into the system: this Q-ball most likely forms from the perturbation first entering 
the lattice, i.e. the first unstable mode which fits inside the lattice. However, the most 
amplified mode is even larger. These considerations make the gauge mediated scenario 
ineffective to study on a lattice: the lattice size should be huge in order to contain the 
relevant modes visible already at the very beginning. 

The situation is, however, different in the gravity mediated case. The instability 
band is characterized by 

^<0'~ V\ (8) 

which has to be positive in order to have growing modes at all. Moreover, the max- 
imally growing mode can be identified as {kmax/iY — ~ ^") growth is 
characterized by l^^l ~ e°''^^\ where 

ait) = ^. (9) 



If ^ ~ m the wave number of the maximally growing mode is k^ax = y ^\K\ a{t) m and 
d; = I |K| m. This would mean that the criteria for the end of linear growth, ^ ~ 1, 
implies a{t) ~ ~ In ^^^^^^^ which, in the light of the simulations, is typically too 
small, as will become evident later. Therefore it is important to recognize that for the 
a; > 1 cases 6{t = 0) <^ m, and hence no instabilities occur because V" is typically 
positive. Due to the dynamics 6 begins to grow and unstable modes appear when 
^2 _ y// changes its sign to positive. At that time, however, the rate of perturbation 
growth is very slow because a is very small [6 ^ m). So the growth of the perturbations 
is very slow at first, but later the rate of growth increases and becomes much faster 
when the average 6"^ is large enough. This lengthens remarkably the time needed to 
enter the non-linear era of perturbation growth and pushes the actual formation of 
Q-balls later. 

Note that when the growth of linear modes is finally speeded up, the most amplified 
modes are remarkably smaller than in the gauge mediated scenario. This naturally 
implies that the size of the largest Q-ball is much larger in the gauge mediated than 
in the gravity mediated scenario. 

In both gauge and gravity mediated cases the unstable perturbations 5$ remain 
linear as long as they are much smaller than the background AD-condensate field $: 

- « 1. (10) 



On the other hand the charge density perturbations 5q 
smaller or equal than the field perturbations, i.e. 



5q 


< 


5$ 


1 







g($ + 5$) — are always 



fill 



This implies definitely that 
perturbations. 



1 indicates the end of the linear growth era of the 



2.2 Equilibrium ensembles 

After fragmentation, the AD lumps are expected to interact vigorously. Gradually the 
field fragments will settle to the state of lowest energy by emitting and exchanging 
smaller fragments. If the interaction is fast enough compared with the expansion 
rate of the universe, a natural expectation is that the final state should consist of an 
equilibrium distribution of Q-balls and anti-Q-balls in a state of maximum entropy. 
We shall now work out the consequences of such an assumption in the case of gravity 
mediated Q-balls, for which the mass Mq is given approximately by 

Mq^tjiQ. (12) 

Q-ball (anti-Q-ball) distributions N^{Q,p) {N^{Q,p)) are subject to the following 
constrains: 

Etot = E+ + E^ , E^= JdQ dp EiQ,p) (Q, p) 

Qtot = Q+-Q-, Q± = IdQdpQN±{Q,p), (13) 

where E{Q,p) ^ A/p^~+m^Q^ is the energy of a single Q-ball, E^ {EJ) and (Q-) 
are the energy and charge of Q-balls (anti-Q-balls), and Etot and Qtot are respectively 
the total energy and charge of Q-balls and anti-Q-balls, which are equal to the cor- 
responding values of the AD-condensate in the beginning unless significant amounts 
of energy and/or charge are transformed into radiation. It then follows from Eq. (0) 
that 

Etot > m (g+ + g_) >mQtot. (14) 

This condition is independent of the precise form of the Q-ball distributions. It is easy 
to see that the energy-to-charge ratio (see Eq. (|)) x > 1. If x = 1, the inequalities in 
Eq. (|T^) simplify to equalities and this gives directly the Q-ball distributions 

N+{Q,p) = N+iQ) 5ip) , N.iQ,p) = , (15) 



which means that there are no anti-Q-balls and that the Q-balls for at rest with 
a distribution N+{Q). However, this is not a realistic case. Indeed, if all the baryon 



asymmetry resides in Q- and anti-Q-balls, then at times earlier than 10~^s, Qtot/Q+ ~ 
AB ~ 10^*^. (Since B — L is conserved in the MSSM, this holds also for the purely 
leptonic flat directions.) From Eq. ( [I^ ) it then would follow that 

X ~ 10^ . (16) 



Even if all of the baryon asymmetry were not carried by Q-balls (or rather, in this 
case, by B-balls), it is likely that some of it is. As a consequence, a natural expectation 
is a; 1. 

Hence the number of Q-balls, A^+, and the number of anti-Q-balls, A^_, are approx- 
imately equal so that the total number of Q-balls is Nfot ~ Then from Eq. (0) 
it follows that 

Ei ^ E{Q,p) JdQdpNi{Q,p) ^ E{Q,p)N+, 

Q± ^ Q j dQ dp N^{Q,p) ^ Q N+, (17) 

where E{Q,p) is the average energy of one Q-ball with the average charge Q and 
average momentum p. Now the energy charge ratio of Q-balls (anti-Q-balls), x+ = 
E^/mQ^ (a;_ = E_/mQ_), determines the average momentum of Q-balls 



p^mQ ^x\-l>0.lmQ , (18) 

when x± > 1.01. The average velocity will then he < v >= p/E{Q,p) > 0.1. If 
Q-balls were non-relativistic, it would cause a restriction x± — 1 < 0.01, which is not 
natural, as will be seen later. Thus the main bulk of Q-balls, except the largest ones, 
may be expected to be relativistic. 

Since most Q-balls are relativistic it is not surprising that the reaction rate turns 
out to be much larger than the Hubble rate. We find 

^ = ^<"^^"^(g:^' ^''^ 

where we have approximated the cross section by the geometrical cross section with 
R the radius of the Q-ball. The cross section actually depends on the relative phase 



difference between the colliding Q-balls [16| as is largest when the phases are aligned 



however, the order of magnitude is that of the geometrical cross section. 

When X ^ 1, the charge and the momenta of Q-balls and anti-Q-balls are almost 
the same. The reaction rate Eq. (|T^ is at its maximum for p = mQ corresponding 
to x± ~ \/2, which is about the same as the equilibrium value. The average charge 
is Q = q+Vi where g+ is the charge density inside a Q-ball and Vi it's volume. Thus 



we can rewrite (|TU|) as F = 3Vip/2mQR = 3p/2mq+R. In the dense-packing limit, 
p ^ 2mg+x+ ~ 2-\/2mg+, and hence 



R 



(20) 



We approximate the radius of Q-balls by i? ~ \K\ ^/^m ^ Eq. and compare 
this to the Hubble rate H ~ mtQ/t to find that T > H when 



t > 



to 



\K\y^ 



(21) 



Since typically \K\ ~ 0.01 — 0.1 [|12|, the reaction rate will be larger than the Hubble 
rate from the very beginning, and a thermal equilibrium can be expected to be reached. 

Guided by these considerations, let us now find the equilibrium distributions for 
Q-balls in the gravity-mediated case in two and three dimensions. The one-particle 
partition function reads 



f d^x d^p f 



Vd 



dQ e 



(22) 



where p is the chemical potential related to the charge of Q-balls. In D dimensions 
we may integrate Eq. (p2D to obtain a distribution function with respect to charge 



D+l oo 

2 



= J dQ Z,{Q) = 2gVDf3\^£^j J dQ |g|^ e^"^ Kn^{l3m\Q\) , (23) 



where Kn{z) is the modified Bessel function. Then we may integrate Eq. ( p3D to obtain 
an expression that can be split into positive and negative charge parts: 



Zi + zt 

2gVD(3m^+\D + 1)\ 
7rfr(f + 2) 



F(Z} + 2, f + 1; f + 2 : 



{f5m T Ai)^+2 



(24) 



where F{a, b;c : x) is the hypergeometric function. For D = 2 and D = 3 these 
expressions simplify to 



Z^ 
Zt 



gV2 2PmTP 



D 



(25) 



p{5(3^m^ - 2/i2) 
2(/32m2 - p^y ' 



+ 



3p^m^ 



2^2 



/i2)2 



arctan 



/ Pm ± p 
(3m =F p 



D 



(26) 



The total partition function is the sum of these two terms: 

2gV2(3m^ 



D = 2 ■ 



D = 3. 



The grand canonical partition function is given, as usual, by 

oo p.N r/N 



N=0 



{z = e^) with the energy, charge and number of Q-balls given respectively as 

d log Zg 



_ d log Zg „ _ d log Zg 
-L^tot — — — 1 Vtot — t: 1 ^^tot 



dp 



(27) 
(28) 



(29) 



(30) 



Using Eq. (|23D in Eq. (^) the total energy of the Q-ball system is obtained 

oo 

Etot = mNtot J dQ fE{Q) , (31) 



where 



fE{Q) 



6- D 



(2vr)- 



D-l 
2 



(/?m|Q|)^ e'^^ Ki^{(3m\Q\) (32) 



is the energy distribution function per Q-ball scaled with the mass m. The average 
energies of Q-balls and anti-Q-balls are then found from Eqs. (|25|,|26|) to be 



Ntot (/3m±/i)2(3/?2m2^3/?m/x + /i2) 



2Ntot - 



± 



Pm =1= n 



D = 2- 



TT/3 L 2(3^m^ 

P2^2 _ 



arctan - 



D = 3. 



/ (5m ± /X 



j3m =F /i 



The total energy is, again, the sum of positive and negative parts: 

iV,,,3/52m2 + /i2 



Etot = 
Etot 



(3 (3'^711? — yU^ ' 

Ntot 4/?'m2 + 



D = 2- 
D = 3. 



(33) 



(34) 

(35) 
(36) 



Similar calculations can also be repeated for charge and particle number. We get the 
total charge 



Qtot = Ntot I dQ sign(g) fgiQ) , 



(37) 



where 



iQiQ) 



2 (27r)- 



D-l 
2 



6-D 



1 - 



{(3m\Q\)^ e^'? Kn+i{Pm\Q\) (38) 



is the charge distribution function per Q-ball. The positive and negative parts of the 
total charge are: 



V± - 1-37-^ o , V - Z ; 



<5± = 



Ntot 



f3m ^ jj, 



± — 1^ t: arctan ' 



/3m ^ II 



(39) 

, D==3 . 
(40) 



The charges are defined to be positive so that the total charge is Qtot — Q+ — Q-- 
Thus 



Qtot — Ntot 
Qtot — Ntot 



4fi 



5/i 



L> = 2; 
D = 3 . 



Moreover, the relative number distribution function of Q-balls fulfilling condition 



(41) 
(42) 



CO 

(3m I dQ fN{Q) = 



(43) 



reads now 



fN{Q) 



2(27r)- 



2 



6-D 



1 - 



(3m 



{Pm\Q\)^ e'^^ Kn+i{Pm\Q\) . (44) 



Finally, the number of positive and negative charge Q-balls is then found to be 
Ntot 



N± 
N± 



AP^m^ 

Ntot 



{2pm^ li){pm±ii)\ D = 2; 



TT 



u(5^W-2^ r—- pm±n 

± tttta — -A \ p^m'^ — + 2 arctan W -; 



, D = 3 



(45) 



(46) 



We have plotted the charge and number distribution functions for — 10 < (3mQ < 
10 in D = 2 and D = 3 for ^ = 0, 0.1(3m in Figs. |l| and |. The energy distribution 
function has a behaviour which is quahtative similar to the number distribution func- 
tion. The total numbers of Q-balls and anti-Q-balls in D = 2 and D = 3 are plotted 
in Fig. ^for 1 < X < 30. When x > 30 there are almost equal numbers of Q-balls and 
anti-Q-balls. 

When total energy and charge are fixed, the chemical potential fi can be found to 
read 



pm (2x - V4x2 - 3) , D = 2; (47) 

= ^ (5x - V25x2 - 16) , D = 3 . (48) 

Using one-particle partition function Eq. (^2]) we can also calculate the average 
velocity of Q-balls to be 

<^>-i^I w 1'"^^ ^il"^ * " ■ 

Vo —00 —00 



from which we get, using Eq. (^3]), the velocity distribution function, which is plotted 
in Fig. I, 

r(^) KuiPm\Q\) 
^""^^ r(f) i f3m\Q\ Ko^i(3m\Q\) ■ 

The average velocity gives again the hypergeometric function and finally we have 

<v>=-^^(3^m^-^i^ , D = 2- (51) 
4pm ^ 



< V >= — P^m^ - , D = 3. (52) 
Snprn * 

Assuming x ^ 1, we see that the chemical potential will almost vanish because 
/i ~ 0{l/x). This in turn implies that the amounts of positive and negative charges, 
Eq. (|39|) , will be much larger than the total charge, Eq. (PD, which is of the order 
of 0{fi). Note that the average velocity is relativistic: ir/A for D = 2 and 8 /Sit for 
D = 3 but the velocity distribution shows that for small Q-balls the velocity is close to 
one whereas for large Q-balls we get non-relativistic velocities. (Thus the qualitative 
behaviour in these two cases is the same, as can be seen in Fig. ^) Moreover, in the 
large x limit (vanishing fi limit) the energy to charge ratio of Q-balls and anti-Q-balls, 
x± = E±/mQ±, is easily calculated. Using Eqs. (^^) and Eqs. (|39| ,^0|) we obtain 
x± ~ 2, regardless of dimension. 

Using the equilibrium distributions one can verify that the reaction rate is larger 
than the Hubble rate so that the assumptions in this Section are self-consistent. The 



real situation is naturally expected to be much more complex, but as we will see in 
the next Section, from numerical simulations one extracts distributions which are very 
close to the thermal ones discussed here. 



3 Numerical simulations 
3.1 Preliminaries 

Because analytical considerations of the previous Section indicate that Q-ball distribu- 
tions in two and three dimensions behave qualitatively in the same way, it is reasonable 
to expect that two dimensional simulations will also shed light on the more realistic 
three dimensional case. The great advantage is of course, that much less CPU time is 
needed for the simulation. Therefore we have simulated Q-ball formation numerically 
on a 2 + 1 dimensional lattice. The largest lattice size used in simulations was 100 x 100 
(simulations on smaller lattices were run to study lattice size effects). We take as the 
initial condition a uniform AD-condensate with an arbitrary phase uo, 

$ = 006*"^* + (5$ (53) 



with uniformly distributed random noise 5$ added to the amplitude and phase. The 
amplitude ratio of the noise and the condensate field, |(5$|/|0o|, is here C(10^^^). 
The phase of the AD-condensate varies in the range uj = 10°, 10"^, 10^^ which 
corresponds to x ~ 1,...,10^. The initial amplitude of the AD-field is set to 0o = 
10^ GeV. This is smaller than the actual condensate amplitude in d=4 fiat direction 
[^] but should be large enough for the simulation to encompass all the qualitative 
features of condensate fragmentation. 

The equation of motion of the AD-condensate is 

$ + 3if$ _ v^'l' + m2$[l + — - + ^log(^)] - gH'l^l + = 0, (54) 

where M is a large mass scale, a is the scale factor of the universe and H is the Hubble 
parameter, H = a/a. 

For the simulations the field was decomposed into real and imaginary parts, (p = 
^(01 + i(t>2)- We also rescale field and space-time according to 

if = — , h = — , r = mt, ^ = mx. (55) 
m m 

The parameter values chosen for the simulations were m = 10^ GeV, K = —0.1 
and A = I (the gH^-teim can be omitted since gH^ <^ m^). The universe is assumed 
to be matter dominated so that a = aot^ and H = ^. The spatial lattice unit 



was = 0.1 and the temporal unit Ar = 0.9A^. The calculations were done in 
a comoving volume so that the physical size of the lattice increases with time. The 
initial time is tq = 100 and a = 0.1 x r^. The simulations were run up to ~ 2 x 10^ 
time steps. 



3.2 Outline of numerical results 

We have studied six different cases, uj = 10°, 10^^ 10^^ (x ~ 1, 10^) on differently 
sized lattices. (Note that even though x may cosmologically be as large as 10* the main 
features of the Q-ball formation process are the same for any x ^ 1. Increasing the 
initial energy in the condensate only increases the relaxation time of the condensate 
and hence CPU time required for the simulation.) In all of the six cases the qualitative 
description of the evolution of the AD-condensate at the beginning of the simulation 
is similar. First the charge density of the condensate decreases uniformly in the box 
due to the expansion of the universe (the charge in the comoving volume is constant 
throughout the simulation). No large fluctuations are visible yet at this time and the 
fluctuation spectrum corresponds to the white noise present in the initial conditions. 
As time progresses a growing mode can be seen to develop. White noise is still present 
but the growing mode soon starts to dominate. This process continues until lumps of 
positive charge develop. 

The further evolution of the AD-condensate depends on the initial energy-to-charge 
ratio X of the condensate and hence on the value of u in the initial conditions. If 
X = 1, i.e. the energy-to-charge ratio is equal to that of a Q-ball, the evolution of the 
AD-condensate continues in a similar fashion. As the universe continues to expand 
the lumps of positive charge slowly develop into Q-balls while their spatial distribution 
effectively freezes in the comoving volume. No negatively charged Q-balls, anti-Q-balls, 
are formed in this case. After going non-linear, the lumps just evolve into Q-balls and 
finally freeze due to the expansion of the universe. This is in fact exactly the case 



studied in |jT9|, where the formation of Q-balls was followed in a three dimensional 
simulation. No anti-Q-balls were however observed since x = 1. 

On the other hand, if u is smaller than one so that x ^ 1, as seems more natural, 
the Q-ball formation process is much more complicated. After the positively charged 
lumps have formed, expanded linearly and then developed non-linearly, the extra en- 
ergy stored in them causes the lumps to fragment as they evolve into Q-balls. In this 
process a large number of negatively charged Q-balls is formed. The total charge in 
the negative and positive Q-balls is approximately equal so that the initial charge in 
the condensate is in fact negligible compared to the amount of charge and anti charge 
created. 

As the universe continues expanding, the Q-ball and anti-Q-ball number density 
distributions freeze as they tend towards the equilibrium distributions discussed in the 



previous Section. 



3.3 Evolution of the AD-condensate 

We have plotted the charge density as a 3D-plot and as a contour plot in the comoving 
volume at different time intervals in the case uj = 10~^ in Figures One should 

note that the z-scale of the 3D-plots and the gray scale coloring of the contour plots 
varies from Figure to Figure^. 

In Fig. 1^ the initial small stochastic perturbations can be seen clearly. The fastest 
growing mode stars to dominate in Fig. ^ while the small perturbations are still 
visible. In Fig. ^ the growing modes are large enough for the initial perturbations to 
be no longer visible (they are obviously still present but due to the increasing scale 
they simply can no longer be seen). The linear growth continues in Fig. |^ until non- 
linear growth begins as depicted in Fig. |^. Rapid non-linear growth progresses and 
in Fig. 1^ the fluctuations of the field are already ~ (9(10^) while up to Fig. H they 
are only at the level ~ O{10^^). In Fig. |10| we can see how the lumps of charge are 
dynamically arranged in string-like features. These filaments are still visible in Fig. 
|TT], where the condensate has further fragmented into lumps while the scale of the 
fluctuations has grown to ~ O{10^^). The filament texture disappears from sight in 
the next two figures, Figs |12| and |13|, while negative charge starts to develop (some 
traces of the filaments can still be seen in the distribution of charge lumps). 

Up to this point the cases x = 1 and a; ^ 1 are qualitatively similar. If a; = 1, no 
negative charge develops and the lumps slowly relax into Q-balls while the universe 
expands which finally freezes the distribution. If, however, a; ^ 1 like in the case 
shown here, new qualitative features become apparent after this point. Negatively 
charged lumps start to become visible in Fig. 0, and in Fig. |15| the ratio of positive 
to negative charge approaches one. The lumps of charge are complicated configurations 
that slowly relax into Q-balls and anti-Q-balls. Finally, in Fig. |16|we can see separate 
Q- and anti-Q-balls. The visible ten or so large Q-balls (anti-Q-balls) are only a very 
small fraction of the total number of Q-balls present in the lattice. 

The different time scales of the fragmentation process can also be seen from Figs ^ 



16| . The period of linear growth continues up to r = 2160 after which rapid non-linear 
growth occurs between r = 2160 — 2520. Negative lumps begin to form at around 
r = 3600 while the relaxation process to Q-balls lasts until r ~ 10^. 



The process of Q-ball formation has been illustrated in more detail in Fig. |T7. 
We are again considering a comoving volume and the scale of the figures is varied 
between different panels. The first panel is taken at r = 5400 and the last panel 



at r = 180000. Here a; = 10 ^ as before. In the first panel two large fluctuations 



^Colour coded versions of the figures, where more details can be seen, can be found at 
www.utu.fi/"' t uomul / qballs / 



that both carry large positive and negative charge densities are visible. Their shapes 
are still far from spherical and at this stage the future evolution of the fluctuations 
is far from obvious. In the following panels it can be seen that the initial complex 
configurations more or less develop into charge-anti charge -pairs that fluctuate and 
move with time. Eventually only one large ball remains from each of the two pairs. 
The final fate of the other partners cannot be determined from the Figures because of 
the presence of a large background. It may be annihilated in the process, or it may 
move away from the larger ball. It is clear, however, that the large balls are formed 
from these complex disturbances through a highly complicated relaxation process. 



3.4 Distributions 

The simulations allow us to determine the charge distribution of Q-balls and anti- 
Q-balls at different uj and time intervals. In practice we look for the distribution 
of the local charge maxima and minima since the charge of a ball is approximately 
proportional to these extremes. This is due to the fact that in the potential studied in 
this work, the Q-ball radius has only a weak dependence on the charge and the profile 
can be well approximated with a gaussian Ansatz [|1^. The cumulative distributions 
are plotted in Fig. ^ for the cases u = 10°, u = 10~^ with a detail of the u = 10^^- 
case enlarged. 

As can be seen from the Fig. |18|, the two cases are quite different from each other 
as expected. In the uj = 10" -case {x ~ 1) there form practically only positive balls, 
whereas in the uj = 10"^ -case an approximately equal number of positive and negative 
balls are created. From the previous description of the Q-ball formation process, these 
differing features are expected. In the detail of the uj = 10~^ -case one can also note 
that there is a small 'bump' at zero charge density. 

We can now compare the distributions obtained from simulations to the equilib- 
rium distributions given in Sect. 2. The cumulative distribution function, 
for a thermal distribution in terms of the probability distribution function, g{Q) = 
g/iQ-/3|Q|^^|Q| _|_ -|^^^ reads as (omitting the fugacity of the distribution) 

Q Q 

F{Q,fi,(3)= J dQ'g{Q')= J dQ'e^'^''^\'^\(3\Q'\ + (56) 



However, as can be seen from Fig. |18|, the distribution is not quite an equilibrium 
one but deviates slightly from the thermal one because of the presence of the 'bump' 
at zero charge density. This feature is, however, a lattice artifact, as we will see. 

Let us account for the 'bump' by introducing a new probability density function 

e<"5-'5lel(/J|Q| + l), |Q|>Q„i. 
= ^ 0. IQI < Q„,.. <^'' 



where Qmin is a parameter to be determined by fitting to the numerical results. With 
this choice of the probability distribution function, the cumulative distribution function 
can be calculated to be 



F{Q, /i, /3, Qmin) = y{mm{Q, -Qmin}) + y{Qmm) - y{m&x{Q, Qmin}), (58) 

where 

r)pf^Q-\Q\P nrt 



We can now fit y{Q) to the numerical data. An example is shown in Fig. ITUI where 



F{Q, fi, P,Qram) is fitted to the u = 10^^-cases at r = 180000. y{Q) fits the data 
generally quite well in all the cases where x ^ 1. If = 10° one obviously needs 
to consider a different type of fitting function since in this case no anti-Q-balls are 
present in the distribution. 

To test whether the inclusion of Qmin to the fitting parameters is justifiable we 
have also fitted F{Q,fi,i3) to the data. In all the cases we found that this gives a 
much poorer fit. In Fig. [2^ we have plotted (ao/a)^/3 at different time intervals for 



the uj = 10~^-case. 

The (ao/a)^-factor is included to account for the expansion of the universe, since in 
two dimensions the total energy density of Q-balls in equilibrium behaves as ~ 1/(3 for 
large f3. On the other hand, the total energy is essentially conserved so that a{t)^/l3{t) 
remains a constant of equilibrium. 

From Fig. pO] it can be seen how j3 quickly drops at around r = 3600 — 5400. 



This corresponds to Figures |14| and |T5|. At this time the negative Q-balls start to 
form and the large initial charge lumps begin to fragment. The distribution then 
quickly thermalizes and is well described by the equilibrium distribution. We have 
run the simulation until r = 180000 but in Fig. |20| we have only plotted a part of the 
whole simulated time interval. However, (ao/a)^/5 stays constant until the end of the 
simulation. 

The presence of the 'bump' in the distribution may be attributed to a number of 
factors. To study whether it is a lattice induced effect we ran a number of simulations 
where the number of lattice points was different while keeping the physical sized of 
the fiducial volume constant. After the fitting procedure we found that Qmin decreases 
with the number of lattice points used. Therefore it is likely that the deviation from 
the simple thermal spectrum is just due to lattice effects. This conclusion is further 
supported by the fact that in the gravity mediated case the charge-energy relation 
is approximately linear, E ^ Q, so that there is no preference for the formation of 
large Q-balls over small Q-ball^. In addition, when one considers the fact that the 
^Note, that Qmin is still a large number 10^ so that the linear charge-energy -relation holds. 



lattice spacing increases with time and that smaller balls have smaller radii it is not 
too surprising to see a small deviation from a purely thermal distribution at small 
charge density. 



4 Conclusions 

We have studied the process of Q-ball formation from the AD-condensate in the early 
universe by numerical simulations in two spatial dimensions. The analytical consider- 
ations presented in Sect. 2 reveal that a typical Q-ball is relativistic. This conclusion 
holds if at least part of the baryon asymmetry of the universe is contained in baryon 
number carrying Q-balls. The reaction rate of Q-ball interactions is then much larger 
than the Hubble rate and hence we may expect that a state of maximum entropy 
is likely to be reached. Assuming a thermal Q-ball distribution we showed that in 
a typical cosmological situation the AD-condensate fragments into practically equal 
numbers of Q-balls and anti-Q-balls so that ~ A^_. The amount of charge con- 
tained in these balls is much larger than the net charge, (5+ ~ Q- ^ Qtot- Estimation 
of the average velocity of a Q-ball confirms that Q-balls are typically relativistic. The 
analytical calculations were carried out in both two and three dimensions and we found 
that the shapes of the number density distributions or velocities are only very weakly 
dependent on dimension. These considerations strongly support the assumption that 
2-d simulations can in fact capture the essential features of the real 3-d case. 

Our numerical simulations were carried out on a 2 + 1 dimensional 100 x 100- 
lattice. The initial condition was chosen to be uniform AD-condensate with small 
perturbations. The initial energy-to-charge ratio x was varied from 1 to ~ 10^. We 
found that regardless of x the qualitative behavior of the AD-condensate is similar 
in the beginning of the simulation. The small perturbations start to grow and due a 
complex dynamical process the AD-condensate fragments into several lumps carrying 
positive charge. These are not yet Q-balls but some excited states |I3[| . 

From this point on the further evolution of the condensate depends on the actual 
value of X. If x = 1, the lumps slowly relax into Q-balls and no anti-Q-balls are present. 
This is what we would expect from analytical studies and this is what has been seen 
in 3-d lattice simulations []TB|, [T^. If, however, x ^ 1, the condensate fragments into a 



large number of Q- and anti-Q-balls, as we have shown here. We also showed that the 
distributions of the Q- and anti-Q-balls obtained from the simulations are thermal, 
which again is in accordance with analytical expectations. 



An interesting feature which is apparent in Figs. 11-13 is that the system of AD 



lumps appears to possess some non-trivial spatial structure. First, the large Q- and 



For very small Q this relation is modified and therefore also the equilibrium distribution for those 
values of Q is not given by Eq. (pq). This does not, however, affect our argument. 



anti-Q-balls can be seen to form in pairs which give rise to a positive spatial correlation 
between them. It should be noted that if small balls are also formed in pairs they are 
likely to have a homogenized spatial distribution due to their large velocities. Second, 
there are long range correlations which show up as filament-like structures in Figs. 



T0| - |12| . The texture visible in the middle phase of the non-linear evolution may well 



survive in the spatial Q-ball distribution also at later times, although it no longer 
is visible. Numerical correlation analysis indeed shows that still at the final stage 



of our simulation (Fig. |T3|) the mean distance between both Q- and anti-Q-balls 
deviates by over four standard deviations from the random uniform distribution. The 
deviation is negative i.e. the mean distance between Q-balls is smaller than expected 
for a random distribution. These features may affect observables such as the cosmic 
microwave background or large scale structure and deserve further study. 

The fact that the Q-ball distribution in the cosmologically realistic case, x 3> 1, 
tends toward a thermal one implies, like the analytical treatment shows, that a large 
number of Q-balls move relativistically. The very large Q-balls that are clearly visible 
in Fig. |16| appear immobile but there is also a large number of smaller fast moving 
Q-balls which can thermalize the distribution. In a thermodynamical sense one can 
view the fragmentation of the AD-condensate into Q-balls as a relaxation process 
where entropy is maximized: from this point of view the appearance of the thermal 
distribution is natural. 

In this work we left the gauge mediated SUSY breaking scenario without the at- 
tention it deserves. This is due to a number of facts that make the simulation of the 
gauge mediated SUSY breaking scenario more cumbersome compared to the gravity 
mediated case. In the fragmentation process the spectra of growing modes play an im- 
portant role since the first mode that goes non-linear determines the size of the fastest 
growing disturbance in the condensate. The SUSY breaking mechanism, i.e. whether 
it is gauge or gravity mediated, determines the size of the fastest growing mode via the 
actual form of the scalar potential. As we have seen, the size of the instabilities in the 
gauge mediated case are larger than in the gravity mediated scenario. Therefore one 
needs a lattice with a much larger physical size. Obviously, one can always increase 
lattice spacing but this happens at the cost of losing resolution. Moreover, the large 
instabilities correspond to large Q-balls forming from the AD-condensate which in the 
gauge mediated scenario have thin walls. To avoid numerical instabilities due to the 
thin walls a high resolution is required. Therefore, in order to avoid finite size effects, 
one would need a much larger number of lattice points than in the gravity mediated 
case. 

In addition to the numerical difficulties the gauge mediated scenario does not ben- 
efit form the analytical arguments presented in Sect. 2. This is due to the different 
energy-charge -relations of the two cases. For large Q-balls in the gauge mediated case 



the energy is related to the charge by ~ instead of the hnear relation E ^ Q 
in the gravity mediated scenario. From one can see that the Q-integral diverges 
since the exp(— + /iQ)-term is dominated by nQ in the large Q limit in the gauge 
mediated scenario. Physically the lack of a proper asymptotic equilibrium state im- 
plies that the system tends towards a state where charge is concentrated in large balls 
(and ultimately in a single Q-ball) as noted also in Ref. |2^. However, one should also 
take into account that the more charge is stored in large balls the more there exists 
leftover energy in the system that must be accounted for. The final state of the system 
would then be strongly dependent on the time at which the distribution freezes due 
to the expansion of the universe and hence on the initial conditions. 

To conclude, our results indicate that in the gravity mediated SUSY breaking case, 
the AD condensate fragments into an almost equal number of Q-balls and anti-Q-balls. 
The evaporation times of the Q-balls, the nature of particles they decay into, and 
possible spatial correlations surviving the fragmentation process and the subsequent 
decays of Q-balls, are obviously very important factors that determine the signatures 
left by a possible Q-ball period in the early history of the universe. 
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Figure 1: Charge distributions of Q-balls for /i = (solid lines) and /j, = 0.1(3m (dashed 
lines) in D = 2 (thin lines) and D = 3 (thick lines). 




Figure 2: Number distributions of Q-balls for // = (sohd hues) and // = 0.1/3m 
(dashed hues) in D = 2 (thin lines) and D — 3 (thick lines) . 




Figure 3: Number densities {N-) of Q-balls (anti-Q-balls), curves above 0.5 (curves 
below 0.5), in D — 2 (solid lines) and D — 3 (dashed lines). 
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Figure 4: Velocity distributions in D = 2 (solid line) and D = 3 (dashed line). 
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Figure 5: r = 90 
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Figure 6: r = 450 
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Figure 7: r = 720 
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Figure 13: r = 3060 
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Figure 14: r = 3600 
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Figure 15: r = 5400 
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Figure 16: r = 180000 




Figure 17: A detail ol the Q-ball formation process 
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Figure 18: Cumulative distributions for (a) a; = 10 ^, (b) a; = 10° and (c) a detail of 
the uj = 10~^-case 
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Figure 19: A typical fit, a; = 10 







50000 1e+05 1.5e+05 

X 

Figure 20: m(3 at different time intervals 



